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Abstract 

We calculate analytically the asymptotic form of quasi-normal modes of perturbations 
of arbitrary spin of a Schwarzschild black hole including first-order corrections. We use the 
Teukolsky equation which applies to both bosonic and fermionic modes. Remarkably, we 
arrive at explicit expressions which coincide with those derived using the Regge- Wheeler 
equation for integer spin. Our zeroth-order expressions agree with the results of WKB 
analysis. In the case of Dirac fermions, our results are in good agreement with numerical 
data. 
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Quasi-normal modes govern the response of a black hole to external perturbations. In 
general, they possess a spectrum of complex frequencies due to the leakage of informa- 
tion into the horizon. Their observation will reveal information about the characteristics 
of the black hole. In asymptotically AdS spaces, they are related to properties of the 
dual conformal field theories on the boundary through the AdS/CFT correspondence. In 
asymptotically flat spaces, apart from observational possibilities, interest in QNMs has 
arisen [1-11] because the asymptotic form of the spectrum was shown to be related to the 
Barbero-Immirzi parameter [12,13] of Loop Quantum Gravity [14-18]. The asymptotic 
form of the spectrum normalized by the Hawking temperature is given by 

^ « -(2n + l)7ri + ln3 (1) 

for scalar and gravitational perturbations. This has been derived numerically [19-23] and 
subsequently confirmed analytically [3,5]. The analytical value of the real part was first 
conjectured by Hod [24] based on the form of the horizon area spectrum proposed by 
Bekenstein and Mukhanov [25] . Its value is intriguing in Loop Quantum Gravity, suggesting 
that the gauge group should be SO{2>) instead of the expected SU{2), as the latter would 
lead to ^lu/Th ~ In 2 asymptotically. 

The asymptotic expression (1) was analytically derived and generalized to arbitrary 
integer spin j [5] 

^ -(2n-|-l)7rz-|-ln(l-|-2cos7rj) (2) 

and a perturbative expansion was established [11]. Extending these analytical methods 
to half-integer spin (such as the Dirac field) is not straightforward. Dirac quasi-normal 
frequencies have been calculated numerically [26-35]. Using the WKB method, it was shown 
that eq. (2) is valid for half-odd-integer spin to leading order {0{n)) [36]. By employing the 
Teukolsky equation, which is applicable to both integer and half-odd-integer perturbations, 
it was shown that eq. (2) is valid for perturbations of general spin [37]. 

Here, we derive the first-order correction to the asymptotic expression (2) thus extend- 
ing the results for integer spin [11]. Remarkably, we arrive at the same expression as a 
function of the spin as in the bosonic case, even though our starting point is different: we 
employ the Teukolsky equation [38] whose solutions can be expressed in terms of confluent 
hypergeometric functions, whereas in the bosonic case discussed in [11], the quasi-normal 
modes were obtained by solving the Regge- Wheeler equation [39] whose solutions are given 
in terms of Bessel functions. In the case of massless Dirac perturbations, our results are in 
good agreement with numerical data [35,40]. 

We are interested in the quasi-normal modes of massless perturbations of the Schwarzschild 
black hole whose metric is 

ds-' = -f{r)dt^ + ^ + r-'dnl , /(r) = l-- (3) 
j{r) r 

where we have chosen units so that the black hole mass, radius of the horizon and Hawking 
temperature are respectively given by 

M = l , ro = l , Th = ^ (4) 

Massless perturbations are described by the Teukolsky equation [38] for arbitrary spin, 
including both integer and half-odd-integer values. This wave equation may be brought 
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into a Schrodinger-like form [37], 



^^+u^-V[r{u)]^^{u) = (5) 

written in terms of the tortoise coordinate 

r* = r + ln(l — r) (6) 

The potential is given by 

for a spin-j field of angular momentum i. 

For quasi-normal modes, we demand the asymptotics 

^{n)^e^'^''* , r*^±oo (8) 

so that the wave is outgoing at infinity (r* oo) and ingoing at the horizon (r^, — oo). 
More precisely, we shall define the boundary conditions in terms of the monodromy along a 
closed contour enclosing the black hole singularity and running counterclockwise, following 
[5]. Such a contour may me deformed so that the monodromy receives contributions only 
from the singular points r = 1 (the horizon) and infinity. Near the horizon and at infinity, 
the wavefunction behaves respectively as 

$ ^ _ iyi{'^+ij/'^) ^ ^ ^ j.i{oJ-ij)^i<^r ^g-j 

The deformed contour runs counterclockwise at infinity and clockwise around r = 1. There- 
fore, the monodromy is 

M = e^^('^-^^/^) (10) 

Next, we deform the monodromy contour so it encircles the black hole singularity. In 
this case, it receives a sole contribution from the rotation by — 27r near r = 0. We shall 
calculate the monodromy following [37] . This will yield the asymptotic form of quasi-normal 
frequencies. We shall then calculate the first-order correction to the asymptotic expression 
and compare with numerical results [35,40]. 

It is convenient to introduce the dimensionless coordinate z = tor^. In the complex z- 
plane, this rotation is of angle — 47r, because for small r, z ~ on account of (6). Expandng 
around the black hole singularity (z = 0), the potential (7) becomes 

.<-) 
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It is interesting to compare with the effective potential one obtains from the Regge- Wheeler 
equation [39] in the case of integer spin [5]. In the latter case, the leading term has a double 
pole at z = 0, but the first-order correction is of the same form [11]. 
The wave equation (5) may be written as 

Ho + -^Hi + . . * = (12) 
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where 



1 _ ^ 

2z 



16^2 



Hi 



A 

z3/2 



(13) 



and we used the expansion (11) of the potential. 

We may solve (12) perturbatively by expanding the wavefunction 

* = ^(0) + J_$(i) + . . . 
The zeroth-order wave equation reads 

7^o^r(o) = 

Its solutions are Whittaker functions. Two independent solutions are 



X = 



3j 



(14) 



(15) 



(16) 



These functions are not defined for even-integer spin. This can be remedied by adopting 
the Whittaker functions Wx^n{—2iz) and W^x^n{—2ie~^'^ z) as a basis instead. The latter 
have more complicated rotation properties (around the origin), so we shall work with the 
basis (16) for arbitrary spin and then analytically continue the results to include the case 
of even-integer spin. 

At large z, the wavefunctions (16) behave respectively as 



*i°^(z) ~ A±{-2iz)^e'' + B±{-2ie'^'z)^e~'^ 



where 



- r(i±2^) _ r(l±2M) , 

^ r{i±i, + xf ' ^ r(i±A.-A)' 

In deriving the asymptotic form of the wavefunctions, we used the identity 



M;,,^(2/) = r(l + 2/x)e- 



1 



r(i + /x-A) 



r(i + /x-A) 



and the asymptotic form of the Whittaker function 



(17) 
(18) 

(19) 

(20) 



Rotating around the z = singularity in the complex z-plane by —An changes the two 
wavefunctions (16) by phases given by 



where we used 

four times. 
Let 



M_A,M(e-'"y) = e-'-^h+i^^MxAv) 



M = 



mil rni2 
m2i m22 



(21) 
(22) 

(23) 
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be the monodromy matrix. Prom the asymptotic form (17) and the rotation property (21), 
we deduce the eigenvalues and corresponding eigenvectors of the monodromy matrix, 

m( ] =e-'-'^( i+) , m( t- =e^--f t- ] (24) 



^ B+ J \ B+ J ' \ B- J \ B. ^ 

They completely determine M. In particular, 

"^^^ = A^B_ - A_B^ (''^ 

The other entries can also be found but are not needed for our purposes. After some algebra, 
we obtain 

mil = e""^^-^^ (1 + 2 cos TTj) (26) 

mil should be equal to the monodromy (10) found using the contour that encircles the 
singularities at the horizon and infinity, 

mil = M (27) 

This determines the spectrum of quasi-normal frequencies. Using eqs. (10) and (26), we 
obtain the asymptotic form of the spectrum, 

^ = -(2n + l)7ri + ln(l + 2cos7rj)+C>(l/Vn) (28) 

where the Hawking temperature is given by (4). This generalizes the result from the Regge- 
Wheeler equation for integer spin [5]. For half-odd-integer spin (e.g., for a Dirac fermion), 
this reduces to 

^ = -(2n + l)7ri + 0(l/v^) , j + ^GN (29) 

showing that asymptotically, the real part of the quasi- normal frequencies vanishes. 

Next, we calculate the first-order correction to the asymptotic expression (28) using the 
method of ref. [11]. It follows from eqs. (12) and (13) that the first-order correction to the 
wavefunction (14) obeys the wave equation 

7^0^(1) = -7^i*(o) (30) 

The corrections to the two basis vectors (16) are given respectively by 

^'£^(z) = a±(z)*f (z) - h^{z)^%\z) (31) 

where we introduced the functions 

dz' 



b±iz) = f^'^^f{z')H,^f{z') (32) 

where W = 2/x is the Wronskian. Using the explicit expressions (16) and (13), we may 
write 

A I'"" dz' 

a±{z) = -;^M,,±^(-2zzOM,,^^(-2iz') 

A dz' 

h±{z) = ^Ma,±^(-2zz')Ma,±^(-2z/) (33) 
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The asymptotic behavior (17) receives a first-order correction. Prom eqs. (14), (31) and the 
identity (19), we deduce that for large z, 



where 



1 



and we defined 



h±A^ , B'=[l + — a± B± 



b±B. 



a± = a±{oo) , b± = 6±(oo) 
The monodromy (25) is corrected to 



mil 



ALB' - A' B' 



Using eq. (35), this may be massaged to 

1 



mil mil 



2i . b+A^B_ + b-A+B+ 



(34) 

(35) 
(36) 

(37) 

(38) 



Then the condition (27) imphes a modification to the asymptotic spectrum (28) including 
first-order corrections, 



^ = -(2n-|-l)7ri+ln(l + 2cos7rj) ^ 



2i ^ b+A_B_ + b-A+B^ 

sm 47r// 



■ + 0{l/n) 
(39) 

This is our main result. It should be compared with its counterpart from the Regge- Wheeler 
equation derived in [11] which is valid for integer spin j, 



Th 



(2n + l)-Ki + ln(l + 2cos7rj 
l + i A sin27rj^2A\^/l j\ /l_i 



47r3/2 sin M 
+0{l/n) , j G N 



4 2 



4 2 



(40) 



where A is given in (11). 

To obtain explicit expressions for the first-order spectrum (39), let us concentrate on 
the case of half odd integer spin. In this case, there is a simplification due to the fact that 



as is evident from (18) for the parameters (16). Eq. (14) simplifies to 



Th 



= -{2n + l)7ri + ^ii±!l 4^ 6_ + 0{l/n) 



From eqs. (33) and (36), we have 



6_ 



A_ 
4fi 



l-OO 

/ dz z-^'\M^,^^{-2iz)f 
Jo 



(41) 



(42) 



(43) 



For half odd integer spin, the Whittaker function in (43) may be written in terms of a 
Laguerre polynomial as 

which simplifies the calculation of the integral (43). However, one ought to be cautious 
with such simplifications, because in general these integrals may only be defined by analytic 
continuation to the desired spin j. 

For j = 1/2 (Dirac fermion), we obtain from (43) 



-(-2i 
J 



Jo 



= 7(-2^)""r(i-0+... (45) 

where the dots represent terms that vanish as j ^ ^. In the limit j = ^, we obtain 

^.^^ = 2A{-2iy/^r Q (46) 
Using eqs. (11), (18), (42) and (46), we deduce the spectrum of a massless Dirac fermion, 



Th 



(2n + l)7ri + ^ (^+^1 r^(^)+0(l/n) (47) 



,=1 

J 2 



To compare with numerical results [35,40], we shall work with the first-order correction 



which is the difference between the first-order (47) and zeroth-order (29) expressions mul- 
tiplied by the Hawking temperature (4). In fig. 1, we compare our analytic expression with 
the numerical results of ref. [35]. We obtain fairly good agreement with the discrepancy 
between numerical and analytic results increasing with angular momentum. This is be- 
cause the first-order correction to the asymptotic form (29) is proportional to the square 
of the angular momentum quantum number i + j. Moreover, it should be noted that the 
corrections are 0{l/y/n), so they become significant for overtones with n < 100. 

In figures 2 and 3 we compare our expression (48) with numerical results for high 
overtones (n > 100) [40] for ^ -|- j = 1, 2, respectively. The agreement improves greatly, as 
expected. In the case of lowest angular momentum {i + j = 1), the agreement is excellent 
(fig- 2). 

For a massless field of spin j = 3/2, we obtain from (43) 



5_ = :^^_2if/2 rdzi-2iz)-"^-h^''(l + -^^ + ..) 
J Jo \ l~j/2 / 



2 



'3 _ l-i 



/ 2(1-j72)2 +••• ^^^^ 
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where we once again omitted terms which vanish as j — > | . Evidently, in the Umit i — ^ | , 



._3=0 (50) 



6_ 

showing that there are no 0{l/y/n) corrections to the spectrum, 

-{2n+l)m + 0{l/n) (51) 

For a massless field of spin j = 5/2, we obtain from (43) 

= '"<-^") ' --'"(^^w72- (l-,-/2)(2-,-/2) +-j 

4 (41 - 40i + 8j2)r f i - f ) 

In the limit j ^ 5/2, we obtain 



(53) 



Using eqs. (18), (42) and (53), we deduce the spectrum of a massless fermion of spin j = 5/2, 



Th 



^ = -{2n + + ^ Q ) + Oil/n) (54) 



"2 



where A is given in (11). 

It is interesting to note that all spectra derived above for fields of half odd integer 
spin may also be obtained by analytically continuing the general expression (40) which 
was obtained for integer spins [11]. All cases we checked agree with (40) but we have been 
unable to find a general argument applicable to arbitrary half odd integer spin. The same is 
true for integer spin spectra derived from the Teukolsky equation, although the calculation 
is more cumbersome (except for zero spin, in which case the Teukolsky equation reduces to 
the Regge- Wheeler equation). 

In conclusion, we have derived analytic expressions for quasi-normal frequencies of 
perturbations of Schwarzschild black holes of arbitrary spin including both bosonic and 
fermionic fields. Our results were in good agreement with numerical data [35,40]. Our 
calculation was based on the monodromy argument of ref. [5] extended to the Teukolsky 
equation [37]. We derived the asymptotic spectrum including first-order corrections ex- 
tending the argument in [11]. We obtained explicit expressions for various spins and found 
that they were in agreement with the analytic expression derived in [11] for bosonic fields 
(after a simple analytic continuation of the spin to non-integer values) . This is remarkable 
given that our starting point was the Teukolsky equation whereas the bosonic spectrum was 
derived in [11] starting from the Regge- Wheeler equation. It would be interesting to further 
investigate the relation between the two equations and understand the range of applicability 
of the Teukolsky equation as the latter offers significant calculational advantages. 



8 



Acknowledgments 

We wish to thank Fu-Wen Shu for useful discussions and R. A. Konoplya for sending us 
his results. G. S. was supported in part by the US Department of Energy under grant DE- 
FG05-91ER40627. The work of S. M. was partially supported by the Thailand Research 
Fund. S. M. also gratefully acknowledges the hospitality of the Department of Physics and 
Astronomy at the University of Tennessee where part of the work was performed. 

References 

[1] O. Dreyer, Phys. Rev. Lett. 90 (2003) 081301; gc-qr/0211076. 

[2] G. Kunstatter, Phys. Rev. Lett. 90 (2003) 161301; gr-qc/0212014. 

[3] L. Motl, Adv. Theor. Math. Phys. 6 (2003) 1135; gr-qc/0212096. 

[4] A. Corichi, Phys. Rev. D67 (2003) 087502; gr-qc/0212126. 

[5] L. Motl and A. Neitzke, Adv. Theor. Math. Phys. 7 (2003) 2; hep-th/0301173. 

[6] E. Berti and K. D. Kokkotas, Phys. Rev. D68 (2003) 044027; hep-th/0303029. 

[7] R. A. Konoplya, Phys. Rev. D68 (2003) 024018; gr-qc/0303052. 

[8] A. Maassen van den Brink, J. Math. Phys. 45 (2004) 327; gr-qc/0303095. 

[9] A. R Polychronakos, Phys. Rev. D69 (2004) 044010; hep-th/0304135. 

[10] E. Berti, V. Cardoso, K. Kokkotas and H. Onozawa, Phys. Rev. D68 (2003) 124018; 
hep-th/0307013. 

[11] S. Musiri and G. Siopsis, Class. Quant. Grav. 20 (2003) L285; hep-th/0308168. 

[12] J. F. Barbero G., Phys. Rev. D51 (1995) 5507; gr-qc/9410014. 

[13] G. Immirzi, Nucl. Phys. Proc. Suppl. 57 (1997) 65; gr-qc/9701052. 

[14] C. Rovelh and L. Smolin, Nucl. Phys. B442 (1995) 593; gr-qc/9411005. 

[15] A. Ashtekar and J. Lewandowski, Class. Quant. Grav. 14 (1997) A55; gr-qc/9602046. 

[16] C. RoveUi, Living Rev. Rel. 1 (1998) 1; gr-qc/9710008. 

[17] T. Thiemann, gr-qc/0110034. 

[18] C. Rovelh and P. Upadhya, gr-qc/9806079. 

[19] S. Chandrasekhar and S. Detweiler, Proc. R. Sac. London, Ser. A 344 (1975) 441. 

[20] E. W. Leaver, Proc. R. Soc. London, Ser. A 402 (1985) 285. 

[21] H. P. Nollert, Phys. Rev. D47 (1993) 5253. 

[22] N. Andersson, Class. Quant. Grav. 10 (1993) L61. 

[23] A. Bachelot and A. Motet-Bachelot, Annales Poincare Phys. Theor. 59 (1993) 3. 

[24] S. Hod, Phys. Rev. Lett. 81 (1998) 4293; gr-qc/9812002. 

[25] J. D. Bekenstein and V. F. Mukhanov, Phys. Lett. B360 (1995) 7; gr-qc/9505012. 

[26] S. Iyer and C. M. Will, Phys. Rev. D35 (1987) 3621. 

[27] V. Cardoso and J. P. S. Lemos, Phys. Rev. D63 (2001) 124015; gr-qc/0101052. 



9 



[28] H. T. Cho, Phys. Rev. D68 (2003) 024003; gr-qc/0303078. 

[29] R. A. Konoplya, Phys. Rev. D68 (2003) 024018; gr-qc/0303052. 

[30] A. Zhidenko, Class. Quant. Grav. 21 (2004) 273; gr-qc/0307012. 

[31] W. Zhou and J.-Y. Zhu, Int. J. Mod. Phys. D13 (2004) 1105; gr-qc/0309071. 

[32] J. Jing, Phys. Rev. D69 (2004) 084009; gr-qc/0312079. 

[33] K. H. C. Castello-Branco, R. A. Konoplya and A. Zhidenko, Phys. Rev. D71 (2005) 
047502; hep-th/0411055. 

[34] H. T. Cho, Class. Quant. Grav. 22 (2005) 775; gr-qc/0411090. 

[35] J. Jing, Phys. Rev. D71 (2005) 124006; gr-qc/0502023. 

[36] H. T. Cho, Phys. Rev. D73 (2006) 024019; gr-qc/0512052. 

[37] 1. B. Khriplovich and G. Yu. Ruban, Int. J. Mod. Phys. D15 (2006) 879; 
gr-qc/0511056. 

[38] S. A. Teukolsky, Phys. Rev. Lett. 29 (1972) 1114; 
S. A. Teukolsky, Astrophys. J. 185 (1973) 635; 
R. Giiven, Phys. Rev. D22 (1980) 2327. 

[39] T. Regge and J. A. Wheeler, Phys. Rev. D108 (1957) 1063. 
F. J. ZeriUi, Phys. Rev. D2 (1970) 2141. 

[40] R. A. Konoplya, private communication. 



10 



numerigal 
analytic 




Figure 1: Quasi-normal frequencies of massless Dirac fermions for various values of angular 
momentum (from top to bottom: i + j = 1,2,3,4,5); solid lines are graphs of our analytic 
expression (48); diamonds represent numerical data [35]. 
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Figure 2: High overtones {n > 100) of massless Dirac fermions for i + j = 1: solid lines are 
graphs of our analytic expression (48); diamonds represent numerical data [40]. 
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Figure 3: High overtones {n > 100) of massless Dirac fermions for i + j = 2: solid lines are 
graphs of our analytic expression (48); diamonds represent numerical data [40]. 
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